Competition between Vortex Unbinding 
and Tunneling in an Optical Lattice 



M. A. Cazalilla,^ A. lucci,^ and T. Giamarchi-^ 

^ Donostia International Physics Center (DIPC), Manuel de Lardizabal 4, 20018-Donostia, Spain. 
^ DPMC-MaNEP, University of Geneva, 24 Quai Ernest Ansermet CH-1211 Geneva 4, Switzerland. 

(Dated: February 6, 2008) 

We study a system of two-dimensional Bose gases trapped in minima of a deep one-dimensional 
optical lattice potential. Increasing the tunneling amplitude between adjacent gases drives a decon- 
finement transition to a phase where coherence is established between neighboring two-dimensional 
gases. We compute the signature of this transition in the interference pattern of the system as well 
as in its rotational response, which provides a direct measurement of the superfluidity in the system. 
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A recent experiment by the ENS group [H has pro- 
vided fairly direct evidence for the mechanism driv- 
ing the remarkable phenomenon known as Bcrczinsky- 
Kosterlitz-Thouless (BKT) transition. The latter was 
predicted in the 1970s and, in an two-dimensional 

(2D) interacting Bose gas, it is caused by the proliferation 
of (thermally excited) unbound pairs of vortices and anti- 
vortices above a temperature, Tbkt, which depends on 
the microscopic parameters of the gas. As the transition 
takes place in a 2D system, long range order is absent [1| 
both above and below Tbkt because of the enhanced ef- 
fect of thermal fluctuations. Thus, the BKT transition, 
albeit continuous, is not characterized by the emergence 
of an order parameter, and therefore does not fit into 
the general theory for phase transitions established by 
Landau in the 1930s. 

In the ENS experiment [H two independent, harmoni- 
cally trapped, 2D Bose gases were created after loading 
a (several tens of uK) cold ®^Rb gas in the standing wave 
formed by two counter- propagating laser beams {i.e. a 
ID optical lattice). Since each 2D gas is a finite-size pan- 
cake, what was observed is more a sharp BKT crossover 
than a phase transition. However, by direct imaging of 
the expanded gas clouds after release from the trap, it 
was possible to show that, at high enough temperature, 
the 2D Bose gases contain free vortices that leave a dis- 
tinctive signature in the interference pattern observed in 
time-of- flight experiments. This provided sufficient direct 
evidence for the BKT mechanism, which could not be 
obtained in other systems, such like liquid helium films, 
where the BKT transition was observed much earlier Q . 

Concerning the possibility of observing this transition 
in trapped ultracold atomic gases, there has been much 
debate in the literature BBlEEIiil. One of the issues 
requiring detailed investigation was that, being the ex- 
perimental systems harmonically trapped, they are inho- 
mogcneous. This, for instance, causes a non-interacting 
2D Boseeas to become a true Bose-Einstein condensate 
(BEG) [3 [l3|, something that is not possible for the 
homogeneous 2D system [1, [l^l. However, recent the- 
oretical calculations Uj have clarified this issue showing 



that, in spite of the inhomogeneity, thermal activation of 
vortex-anti-vortex pairs, and thus a BKT-like crossover, 
can take place also in the trapped systems. This was 
later confirmed by the experiments [l|. 

As we show in this paper, ID optical lattices, like the 
one used in 1] , offer the possibility of exploring a much 
richer kind of phenomena, which lies beyond the stan- 
dard BKT physics: due to the large degree of tunabil- 
ity of these atomic systems, it is experimentally possi- 
ble to control the tunnelling amplitude between the pan- 
cakes that are formed in the deep wells of the ID lattice. 
This adds a new dimension to the phase diagram of the 
system: the Josephson coupling between the 2D gases, 
which indeed can be more easily controlled than temper- 
ature. When the amplitude of this coupling, is made 
large enough so that tunneling of an appreciable num- 
ber of bosons takes place over the typical duration, tcxp, 
of an experiment {i.e. h/toxp t±), but small enough 
so that the pancakes remain 2D {i.e. tj^ <C /i, /i being 
the chemical potential of the 2D gas), we find that the 
Josephson coupling can drive a transition where phase 
coherence builds up between neighboring pancakes. The 
competition between the interplane phase coherence and 
the thermally excited vortices in the pancakes leads to 
a transition analogous to the deconfinement transition 
studied by some of us [13, [13, [l3| in various types of 
coupled quantum ID interacting Bose systems, such like 
those that are created in anisotropic optical lattices [13] ■ 
Universality between classical statistical mechanics in 2D 
and quantum field theories in ID allows to use the same 
techniques [13, [l^, [l6|, (l3| to tackle the coupled classical 
two dimensional systems as well. In addition to their in- 
trinsic interest, the physics of coupled 2D Bose systems 
is also relevant to the theoretical understanding of re- 
cent experiments in the high- Tc cuprate superconductors 
(see [16| and references therein), for which control over 
the system parameters is much more limited. 

As shown in what follows, when only two pancakes 
are coupled, the resulting system behaves, below a cer- 
tain temperature Tc (see lower panel of Fig. [1]), as sin- 
gle 2D superfluid. Furthermore, if an infinite number 
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of pancakes are coupled, the deconfinement transition 
yields a true BEG. Here we focus on the experimental 
consequences of this transition: although in the BEG or 
superfluid phases the system displays true interference 
fringes, a measurement of the interference contrast like 
the one conducted by the ENS group will yield a jump, 
Aa = i, in the exponent (a) that characterizes the scal- 
ing behavior of the interference contrast with the imaged 
area (Aa = i for the BKT transition [E Ell). How- 
ever, if the interference pattern serves as probe of the 
coherence between the pancakes, it is not a probe of the 
superfluid properties of phase-coherent regime. We show 
that a direct probe of the superfluidity is provided by 
measuring the moment of inertia of the system at low 
rotation frequencies. In the presence of weak Josephson 
coupling between the pancakes, the rotational response 
exhibits a behavior that differs from that of a single {i.e. 
uncoupled) pancake (see upper panel of Fig.[T]). 

In the regime where phase fluctuations are classical, 
i.e. for T > fi, the Bose field in the m-th (m = 1, . . . , N) 
pancake of a ID lattice can be written (away from the 

vortex cores) as "^mir) = Po^^ e*®'"^'') (r = {x,y)). Thus 
the ID lattice is described by the following Hamiltonian 
for a stack of N coupled XY models: 
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where i?xY = ^J2m=i I dr (yQmi'''))'^ is the Hamil- 
tonian of the uncoupled pancakes of area A. We have 
assumed the pancakes to be identical and, for most of 
what follows, we shall take them to be infinite {i.e. 
A — » oo). The bare parameters are K{0) — Trffpo/MT, 
where M is the boson mass, T the temperature, and 
po = Po{T) — |^'m(r)p, the quasi-condensate den- 
sity d, Q ; the (dimensionless) strength of the Josephson 
coupling gj{0) — TTaQPQt±/T, where ag is of the order 
of the healing length ^ = yjh? /Mp. It can be shown 
that, close to the BKT transition, the partition function 
of ((1]) can be obtained from the partition function of a 
ID quantum model with the following Hamiltonian [l9| : 
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The operators 9^ and dx(f>m are canonically conjugated 
fields. The coupling g„ = 27re~^'^/-^, where Eq — 
1.6n fPpo{T)/M is the vortex core energy [Uli^l. 

To ascertain the nature of the phases of the ID op- 
tical lattice, we apply the renormalization group (RG) 



to ([2]), i.e. we iteratively coarse-grain short-distance de- 
grees of freedom to obtain a simpler description of the 
system. To analyze the dependence on the number of 
coupled pancakes, we consider two limiting cases: the 
double pancake system (A'' = 2) and the infinite lattice 
{N ~ oo). We show that the main features of the transi- 
tion in both cases are similar. In the calculations, we have 
exploited that ^ also describes a system of Josephson- 
coupled quantum ID Bose gases in a periodic potential 



at commensurate filling 
RG procedure can be found in [1 

Let us start with the two pancake system 




._More details on the 
We intro- 



duce symmetric and anti-symmetric field combinations 
[13 (/)± = ((/-i ± ^2) /V2 and 0± = {01 ± 62) /V2. The 
action for these fields depends on the parameters K± , 
and gj, which flow under the RG transformations as de- 
scribed by the following set of equations 
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The above equations describe the crossover between two 
strong coupling RG fixed points: when the bare couplings 
determined by T,pQ,t±, and M, are such that g„ be- 
comes of order one first (and the perturbative method 
employed to obtain (|4|) breaks down), vortex unbinding 
takes place before phase coherence between the two pan- 
cakes is established. In this phase the two pancakes are 
decoupled and behave as normal Bose gases. On the 
other hand, when the bare couplings are such that gj 
becomes of order one first, the Josephson coupling term 
(cx gj) causes the relative phase 9- of the two pancakes to 
lock to the same value (modulo 2tt). In this (superfluid) 
phase, only the symmetric field 9+ can fluctuate at large 
distances. This causes the same decay of correlations in 
both pancakes: (V'm(r-)'0m(O)) ~ \r\-^/'^^+ (to = 1,2). 
We note that is related to the superfluid density, 
Ps{T). The relationship is obtained from the transforma- 
tion of ^ under a phase twist where dx9m dx9m + ol 
{i.e. dx9+ dx9+ + ay/2 but 9^ 9^). Thus, 
K+/TT = h^p,{T)/MT. 

The case of large number of pancakes N 00 can be 
treated by a similar application of the RG to dp. This 
leads the following set of flow equations 13, M 113] : 
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^ = {2-K)g^,^=(2-±-j^lgj. 
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The physics described by these equations is similar: as 
we coarse grain the short distance degrees of freedom of 
the system, gj and g^j get renormalized. For given val- 
ues of T, t±, Po, and M, whichever gy or gj reaches a 
value of order one first, determines the phase of the sys- 
tem. When the coupling gy dominates, vortex unbinding 
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destroys phase coherence within the pancakes and there- 
fore suppresses Josephson tunneUing at long distances: 
the pancakes decouple and an array of normal gases is 
obtained. On the other hand, if gj dominates over g^^ 
phase coherence builds up and a BEG forms. The result- 
ing system is an anisotropic 3D BEG, whose properties 
can be computed by using Landau-Ginzburg theory 14 1 . 
By the same 'phase-twist' argument given above, it is 
found that KJtt = h^p,(T)/MT. 

Using the above RG equations we can compute the 
observable consequences of the deconfinement transition. 
Let us first discuss the interference contrast that was 
measured by the ENS group [J, It was argued 

in Ref. [18.] that interference between two independent 
identical Bose systems provides a way to probe corre- 
lations in each system. This idea was employed by the 
ENS group to observe the BKT transition by creating a 
double-pancake system in a ID optical lattice and mak- 
ing the atom hopping between them negligible over the 
time scale of the duration of the experiment. However, 
as discussed above, when hopping between pancakes is 
increased, it leads to a phase where they are not inde- 
pendent anymore. Thus, the expansion images taken af- 
ter release of the atom clouds from the traps will exhibit 
true interference fringes. For instance, for two pancakes 
in the superfluid phase, the imaged (integrated) density 
distribution of the expanded coulds will exhibit oscilla- 
tions that do not wash out after averaging over many 
experimental runs. On the other hand, in the normal 
gas phase (as well as in the superfluid phase at i_L = 0) 
where the pancakes are independent, these oscillations 
average to zero. However, if while seeking for signatures 
of the transition between these two phases, one insists 
in measuring the interference contrast, we show below 
that the latter also provides information about the phase 
correlations between nei ghb oring 2D gases. 

For N = 2, following [18|, we consider the integrated 
interference contrast. 
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where Fi2(ri,r2) = (^I{ri)^^{ri)^\{r2)-^2ir2)); 
(|CqP) is the square of the amplitude of the oscillations 
in the interference pattern with wave number Q = md/t, 
normalized to (the square of) the imaging length {Ly 
is kept fixed) P, [l^ . In terms of symmetric and antisym- 



metric fields, Fi2(ri,r2) 
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In 



the normal Bose gas phase, the relative phase 6- under- 
goes wild fluctuations because vortex-anti-vortex pairs 
unbind. Thus correlations of 0_ asymptotically decay to 
zero exponentially. This implies that (|CqP) ^ L^^" 
with a = ^ 18|. On the other hand, when the Joseph- 



son coupling dominates, long-distance behavior of the 
system, the energy of the system is minimized by lock- 
ing 9- into one of the minima of the cosine in the sec- 
ond term of 0. Hence, in the |ri — r2| oo limit. 
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FIG. 1: Upper panel: Ratio of the moment of inertia to the 
classical moment of interia: I /Id = (1 — ps/p) for N = 2 and 
N ^ oo; t± — 0.05 i-i. To mimic finite size efi'ects we have as- 
sumed that the ID lattice is made of circular pancake-shaped 
2D Bose gases of radius R — 100 5, where ^ is the healing 
length, and therefore we have stopped the renormalization- 
group flow at the scale of R (furthermore, we have assumed 
that po/p = 0.95). Lower panel: Phase diagram (in the ther- 
modynamic limit). We have indicated with an arrow the jump 
in the exponent (a) that describes the scaling of the interfer- 
ence contrast with the imaged length (see Eqs. (|6l8 pV 
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= const. 

Therefore, in the superfluid phase a = 0. Thus, in go- 
ing from the superfluid phase to a phase of decoupled 2D 
normal Bose gases by changing T (or, easier, t±/p), we 
predict a stronger change in the scaling behavior of the 
interference contrast: (|CqP) ~ L~'^°', where a changes 
from in the superfluid to ^ for the normal gas phase. 
Notice that this is a much larger jump than the one pre- 
dicted 18 1 (and observed [l[) for the standard BKT tran- 
sition in a 2D Bose gas, where a jumps from a — ^ at 
T < TxT to a = 



i at T > Ti 



KT- 



For N = oo, the interference pattern receives contribu- 
tions from an infinite number pancakes. Thus, the corre- 
lation function of the integrated density after expansion 
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reads: 



p 

where z is the coordinate of along the axis normal to the 
pancakes and the p-th contrast is 
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(vl'j„+p(ri)vl/„(ri)vl't^(r2)*„+p(r2)). (8) 



Nevertheless, in the (superfluid) BEC phase, the above 
correlation functions also decay to a constant asymptot- 
ically for |ri - ral oo. Therefore, (|C^P) ^ Lf" with 
a = 0. In the normal gas phase, the pancakes decouple, 
and one finds (|Cgp) ~ Li^f with a = \- Thus the same 
jump in a as for N — 2 should be observed. 

Whereas the interference contrast provides evidence 
for the existence or lack of phase coherence between 
pancakes, it is not a direct measurement of the super- 
fiuid density ps . Such a measurement can obtained from 
the response of the system to a slow rotation about the 
axis of the ID lattice (z-axis). In the (rotating) refer- 
ence frame where the gas is stationary, the equilibrium 
state is found by minimizing the free energy F{Q) — 
-T lnZ[r2], Z[n] = Tre-(^-^^^-)/^ with H given by 
(H]) and = J^i J2ap ^apTaiPpi being the z component 
of the total angular momentum operator {ri,pi are the 
position and canonical momentum operators for the i- 
th particle, respectively), and e^y = —^yx = +1- By 
treating the term flL^ as a perturbation (assuming f2 
to be small), we can write ilL^ = J dr A^{r) ■ j{r), 

Jir) = IELi [H.irmrnir)-V^lir)^„^{r)] be- 
ing the momentum density and A^^{r) = fl{—y,x) an 
effective transverse {i.e. V- A^^(r) — 0) vector potential. 
Within linear response theory, the angular momentum 
of the lattice {L^)n=^o = J^ap I dreapVa {jpir))n^o 
Ec/sJ drdr' eo.0rpCo.0{r,r')A^{r') + 0{n^). Thus 
the response to rotation is given by the momentum 
density correlation function at zero rotation frequency, 
Ccf3{r,r') = ^{ja{r)j0{r'))n=n- The latter can be 
decomposed into a longitudinal and a transverse part, 
Cafiiq) = F\\{q)qaqp + F±{q) {5ai3 - QaQis), in Fourier 
space, where qa = Qa/q, and -^11(9) and F±{q) are func- 
tions of g = \q\ only. Since q ■ A (r) = 0, only the trans- 
verse part of the correlation function, which is associated 
with the irrotational part of the gas flow caused by the 
vortices [2l|, will contribute to the response to rotation. 



For small ft, the vector potential is a sufficiently slow per- 
turbation so that only the q = component of the cor- 
relation function matters, and it can be shown [2l| that 
F±{0) = p„ = p — Ps, where p is the gas density. Hence, 
{Lz)n^o = — {l~ Ps/ p)Ici^, where Id is the classical 
moment of inertia of the pancake. Since ps is related to 
(iV = 2) or Ks (N — > 00), the behavior of the ratio 



I /Id = (1 — Ps/p) can be obtained from the RG equa- 
tions (I4l5p . For = 2, in the normal gas phase, 
{Ks 0, for N = 00) and I /Id 1- On the other hand, 
in the superfluid phase, {Ks) flows to a finite value 
and therefore I /Id < 1 because the superfluid fraction 
does not respond to rotation. The transition between 
these two phases is shown in the upper panel of Fig. [1] 
The (abrupt, for = 1) change in I/Ic between the su- 
perfluid (or BEC) and the normal gas phases is rendered 
smooth(er) by finite-size effects, which are mimicked by 
stopping the RG flow at a length scale equal to the size 
of the pancake {R = 100^ in Fig. [J). The change may be 
further smoothen by inhomogeneity effects [lO|. Exper- 
imental determination of I /Id by measuring (i2)n^o is 
not an easy task. Nevertheless, obtaining information on 
{Lz)n, and hence on I / Id, is possible by measuring the 
shape oscillations of the rotating pancakes, the so-called 
'scisors mode' oscillation caused by a sudden rotation of 
a quadrupolc deformation of the trap [2^ . 
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